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NEUTRINO MASS AND MIXING PARAMETERS: A SHORT REVIEW 



G.L. Fogli, E. Lisi a , A. Marrone, A. Palazzo, A.M. Rotunno 
Dipartimento di Fisica and Sezione INFN di Bari, 
Via Amendola 173, 70126 Bari, Italy 



We present a brief review of the current status of neutrino mass and mixing parameters, 
based on a comprehensive phenomenological analysis of neutrino oscillation and non-oscillation 
searches, within the standard three-neutrino mixing framework. 
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1 Introduction 

There is compelling experimental evidence^ that the three known neutrino states with definite 
flavor v a (a = e,/U,r) are linear combinations of states with definite mass (i = 1,2,3), and 
that the Hamiltonian of neutrino propagation in vacuum^ and matter^ does not commute with 
flavor. The evidence for flavor nonconservation (i.e., "neutrino oscillations") comes from a series 
of experiments performed during about four decades of research with very different neutrino 
beams and detection tec hniq ues: the solar neutrino ^ exp eriments Homestake^, Kamiokande^, 
SAGE P GALL EX-GNqEEI, Super-Kamiokande (SK)EI an d the Sudbur y Neut rino Observatory 
the long-baseline reactor neutrin o experiment 

KamLANrjECSl 

the atmospheric 

neutrino experiments Super-Kamiokande ^^ ^ , MACRO ^^j, and . Soud an- 2 and the long- 
baseline accelerator neutrino experiment KEK-to-Kamioka (K2K) ^^l 

Together with the null results from the CHOOZ^l short-baseline reactor experiment, the 
above oscillation data provide stringent constraints on the neutrino mixing matrix, on the split- 
tings between squared neutrino masses, and on matter effects. The absolute neutrino masses 
are being probed by different, non-oscillation searches: beta decay experiments'™, neutrinoless 
double beta decay searches (0^2/3) E^ESl an( ^ precision cosmology ^^^^^^^^^^ i Current non- 
oscillation data provide only upper limits on neutrino masses, with the only exception of the 
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Heidelberg-Moscow 01/2(3 experiment!^, whose claimed signal implies a lower bound on neutrino 
masses. 

Basically all the data are consistent with the simplest extension of the standard electroweak 
model needed to accommodate nonzero neutrino masses and mixings, namely, with a scenario 
where the three known flavor states i/ C(1T are mixed with only three mass states ^1,2,3, no other 
states or new neutrino interac tion s being needed. This "standard three-neutrino framework" 
(as recently reviewed, e.g., in 1221, where further references can be found) appears thus as a 
new paradigm of particle and astroparticle physics, which will be tested, refined, and possibly 
challenged by a serie s of new, more sensitive experiments planned for the next few years or even 
for the next decades ^ j. T he first challenge might actually come very soon from the running 
MiniBooNE experiment!^, which is probing the only piece of data at variance with the standard 
three-neutrino framework, namely, the controversial result of the Liquid Scintillator Neutrino 
Experiment (LSND)^. 

In this short review we focus on the current status of the standard three-neutrino framework 
and on the neutrino mass and mixing parameters which characterize it. The parameters are 
defined as follows. The unitary mixing matrix U, in terms of one-particle neutrino states |i>), is 
defined as (see, e.g.,^), 



K> = E^k>- (i) 

i=l 

A common parameterization for the matrix U is: 

u = o 23 r s o 13 rlo 12 , (2) 

where the O^'s are real Euler rotations with angles 8ij G [0, vr/2], while r,5 embeds a CP- violating 
phase S 6 [0, 2n], 

r 5 = diag(l,l,e +i5 ) . (3) 

By considering r^O^r! as a single (complex) rotation, this parametrization coincides with the 
one recommended [together with Eq. (fTjl] in the Review of Particle Properties^. 

For the sake of simplicity, the phase 5 will not be considered in full generality hereafter. 
Numerical examples will refer only to the two inequivalent CP-conserving cases, namely, e lS = 
±1. In these two cases, the mixing matrix takes a real form Uqp, 



1 C13C12 S12C13 ±«13 

-S12C23 T S23S13C12 C23C12 T S23S13S12 S23C13 j • (4) 
\ S23S12 T S13C23C12 -S23C12 T S13S12C23 C23C13 



where the upper (lower) sign refers to 5 = (S = w). The two cases are formally related by the 
replacement S13 — > — S13. 

The three-neutrino mass spectrum ^^=1^,3 is formed by a "doublet" of relatively close 
states and by a third "lone" neutrino state, which may be either heavier than the doublet 
("normal hierarchy," NH) or lighter ("inverted hierarchy," IH). Here, the lightest (heaviest) 
neutrino in the doublet is called u\ (^2), so that their squared mass difference is 

5m 2 = m| — mf > (5) 

by convention. The lone state is then labeled as ^3, and the physical sign of m| — m\ 2 distin- 
guishes NH from IH. 

Concerning the second independent squared mass difference Am 2 , we define it as ' 



so that the two hierarchies (NH and IH) are simply related by the transformation +Am 2 — > 
—Am 2 . The largest and next-to-largest squared mass gaps are given Am 2 ± Sm 2 /2 in both 
cases. More precisely, the squared mass matrix M 2 reads, in such convention, 



M 2 = diag(m 2 , ml m§) = m 2 + m i i + diag ^- 



5m 2 5m 



-\ — —,±Am 2 



(7) 



where the upper (lower) sign refers to normal (inverted) hierarchy. 

In the previous equation, the term proportional to the unit matrix 1 is irrelevant in neutrino 
oscillations, while it matters in observables sensitive to the absolute neutrino mass scale, such 
as in /?-decay and precision cosmology. In particular, /3-decay experiments are sensitive to the 
so-called effective electron neutrino mass mp, 



|2 2 

rnf 



22 2,22 2,2 2 
c 13 c 12 m l ~r c 13 s 12 m 2 + s 13 m 3 



(8) 



as far as the single i>\ mass states are not experimentally resolvable. On the other hand, precision 
cosmology is sensitive, to a good approximation (up to small hierarchy-dependent effects which 
may become important in next-generation precision measurements'^) to the sum of neutrino 
masses E^^OI 

S = mi + 77i2 + rn?, • (9) 

Finally, if neutrinos are indistinguishable from their antiparticles (i.e., if they are Majorana 
rather than Dirac neutrinos), the mixing matrix U acquires a (diagonal) extra factor 



U^U-Um 



(10) 



containing Majorana phases fa, which are irrelevant in oscillations but not in neutrinoless double 
beta decay (0u2/3). Using the parametrization 

U M = diag (l,e^ 2 ,e^+25)) f 
the expression of the effective Majorana mass mm probed in 0v2[3 experiments^ takes the form: 



U ei m < 



c i3 c i2 m i + c i3 s i2 m 2e + •s 13 m 3 e 



(12) 



Finally, we remark that the constraints on the neutrino oscillation parameters shown here- 
after have been obtained by fitting accurate theoretical predictions to a large set of experimental 
data, through either least-square or maximum-likelihood methods. In both cases, parameter es- 
timations reduce to finding the minimum of a x 2 function and to tracing iso-A^ 2 contours 



around it. We adopt the convention used in 
parameter space obeying the inequality 



1 



and call "region allowed at nc" the subset of the 



Ax 2 < n 2 



(13) 



The projection of such allowed region onto each single parameter provides the na bound on such 
parameter. In particular, we shall also directly use the relation \J A\ 2 = n to derive allowed 
parameter r ange s at n standard deviations. Numerical results and figures are taken from the 
rece nt revie w"^ (to which the reader is referred for details), which also makes use of results from 
Re f s . EHEHI F or previous reviews on neutrino parameters see, e.g., R e f S- l4m 41 | 42 | 43 | 44 | 45 | 46 |4JJ 



2 Constraints on (5m 2 , sin 2 #12, sin 2 #13) from solar+KamLAND data 



It is well known that the angle #13 is relatively small and possibly zero. For #13 = 0, both solar 
and (long-baseline) reactor neutrino oscillations depend solely on the parameters (5m 2 , #12). 
Figure 1 shows the current con straints on such para met ers from a global analysis of all the 
available solar neutrino data^l and of KamLAND data^Sl, both separately and in combination. 
Although (at 3a) multiple solutions can explain KamLAND data, the combination with solar 
data provides a well-defined and unique solution at large mixing angle (LMA) in the mass-mixing 
parameter space. The identification of such solution represents one of the most impressive recent 
advances in neutrino physics. 

Further progress can be expected in narrowing the parameter space in Fig. 1. The 5m 2 
uncertainty is currently dominated by the Kam LAND observation of half-period of oscillations^^ 
and can be improved with higher statistics'-^. The sin 2 9\2 uncertainty is instead dominated by 
the SNO ratio of charged-to-neutral current (CC/NC) event rates, which can also be improved 
with future data^l 

The cur rent solar LMA solution, as compared with results prior to the complete SNO-II 
data set^^, is slightly shifted toward larger values of sin 2 9 12 an d allows higher values of 5m 2 . 
[Our current best-fit point for solar data only is at 5m 2 = 6.3 x 10~ 5 eV 2 and sin 2 9\2 = 0.314.] 
This trend is substantially due to the larger value of the CC/NC ratio measured in the complete 
SNO II phase (0.34™31) 

with respect to the previous central value (0.31^2). We also find that 
the SNO-II charged-current spectral data ^ contribute to allow slightly higher values of 5m 2 
with respect to older results. 

For #13 > 0, the solar and KamLAND v parameter space is spanned by (5m 2 , sin 2 9\2, sin 2 #13). 
Figure 2 shows the current 2a bounds in such space (both separately and in combination) in 
each of the three coordinate planes. Remarkably, both solar and KamLAND data are consistent 
with #13 being small (~ at best-fit), in agreement with independent atmospheric, accelerator 
and short-baseline reactor data (see the next section). The combined upper bound on sin 2 #13 
in Fig. 2 is at the interesting level of ~ 5%. 



3 Constraints on (Am 2 , sin 2 9 23 , sin 2 #13) from SK A tm+K2K-|-CHOOZ data 

In the limit 5m 2 / Am 2 <C 1 (one-dominant-mass-scale approximation), the leading parameters 
in atmospheric and long-baseline accelerator searches are (Am 2 , sin 2 #23 ; sin 2 #13) . Sublea ding 
effects induced by 5m 2 7^ (i.e., L MA effects in terrestrial neutrino 

oscillations™^) 

arc 

present, however, even for #13 7^ 0^"} In accurate calculations, it is worthwhile to include 
such effects numerically, e.g., by fixing (5m 2 , sin 2 #12) at their best-fit value in Fig. 1 in a full 
three-flavor analysis of atmospheric and K2K data, as done in the following two figures. 

Figure 3 shows, for #13 = 0, the results of our analysis of SK (atmospheric) and K2K data, 
both separately and in combination. The K2K constraints are octant-symmetric and relatively 
weak in sin 2 #23 > while they contribute appreciably to reduce the overall Am 2 uncertainty. The 
SK atmospheric neutrino contraints are instead strong on both mass and mixing parameters, and 
also slightly asymmetrical^^ in sin 2 #23- Unfortunately, current data are not accurate enough to 
promote this slight asymmetry to a real #23-octant discrimination. However, it is not excluded 
that future, high-statistics atmospheric neutrino data might be able to do so, if #23 is not too 
close to 7r/4l"*l. Such po ssible #23 octant asymmetry, together with a measurement of #13, is 
crucial for model buildingESEH. 

For #13 > 0, SK+K2K data are also sensitive, in principle, to the neutrino mass hierarchy 
[sign(±Am 2 ) = ±1] and to the CP parity [cos<5 = ±1]. However, the dependence is very small 
within the CHOOZ bounds on #13 (see, e.g., Ref.l^and references therein), and thus it makes 
sense to marginalize the SK+K2K+CHOOZ \ 2 function with respect to hierarchy and CP parity. 



The results are shown in Fig. 4, in terms of the projections of the (Am 2 , sin 2 #23 > sin 2 #13) region 
allowed at 1, 2, and 3a onto each of the coordinate planes (with LMA effects included). The 
best fit is reached for nonzero #13 (mainly due to a slight preference of low-energy atmospheric 
data for v e event appearance), but #13 = is allowed within less than la. The preferred value 
of sin 2 623 remains slightly below maximal mixing. The best-fit value of Am 2 is 2.4 x 10 -3 eV 2 . 
Notice that the correlations among the three parameters in Fig. 4 are very weak. 

4 Global constraints on oscillation parameters 

The results of the global analysis of solar+KamLAND data (Sec. 2) and of SK+K2K+CHOOZ 
data (Sec. 3) can now be merged to provide our best estimates of the five neutrino oscillation 
parameters (5m 2 , Am 2 , 612, O13, O23), marginalized over the 2x2 cases with different mass hierar- 
chies and CP parities (which are physically different but phenomenologically indistinguishable at 
present). The bounds will be directly shown in terms of the "number of sigmas", corresponding 
to the function (Ax 2 ) 1 / 2 for each parameter. 

Figure 5 shows our global bounds on sin 2 #13, as coming from all data (solid line) and from the 
following partial data sets: KamLAND (dotted), solar (dot-dashed), solar+KamLAND (short- 
dashed) and SK+K2K+CHOOZ (long-dashed). Only the latter set, as observed before, gives 
a weak indication for nonzero #13. Interestingly, solar+KamLAND data are now sufficiently 
accurate to provide bounds which are not much weaker than the dominant SK+K2K+CHOOZ 
ones, also because the latter slightly prefer #13 > as best fit, while the former do not. 

Figure 6 shows our global bounds on the four mass-mixing parameters which present both 
upper and lower limits with high statistical significance. Notice that the accuracy of the param- 
eter estimate is already good enough to lead to almost "linear" errors, especially for Sm 2 and 
sin 2 812 ■ For Am 2 and sin 2 623 , such "linearity" is somewhat worse in the region close to the 
best fit (say, within Ha), and thus 2a (or 3a) errors should be taken as reference. 

5 Non-oscillation data and their interplay with oscillation constraints 

Since oscillation data fix the mass splittings dm 2 and Am 2 , the observables sensitive to absolute 
neutrino masses (mp, mgg, X) are high/ correlated with each other, both in normal and in inverted 
hierarchy: typically, when one increases the other also increases. Therefore, upper bounds on 
any of them translate into upper bounds on the others. However, the upper bounds on mp are 
currently weak (a few eV)^l and the relevant discussion can be limited to (mas, S) at present. 

Figure 7 shows the impact of all the available non-oscillation data, taken at face value, in 
the parameter space {m gg, £), at the 2a level. The horizontal band is allowed by the positive 
0v2f3 experimental claim!^ equipped with the nuclear uncertainties oJ^as described id^. The 
slanted bands (for normal and inverted hierarchy) are allowed by all other neutrino data, i.e., by 
the combination of neutrino oscillation constraints (from Figs. 5 and 6) and of astrophysical and 
cosmological constraints fro m C osmic Microwave Background (CMB) , large scale str uctu res 
from galaxy su rveys (2dF)^3, and small scale structures from Lyman a forest data ^| as 
described in The tight cosmological upper bound on S prevents the overlap between the 
slanted and horizontal bands at 2a, indicating that no global combination of oscillation and 
non-oscillation data is p ossi ble in the sub-eV range. The "discrepancy" is now even stronger 
than it was found in Ref.^^, due to the adoption of smaller 0^2/3 nuclear uncertainties'^. It is 
premature, however, to derive any definite conclusion as to which piece of the data (or of the 3v 
scenario) is "wrong" in this conflicting picture. E.g ., b y relaxing either the 0v2(3 lower bound 
or the Lya data, global combinations are possible 1 ^. Further experimental and theoretical 
research is needed to clarify the interplay of absolute neutrino observables in the sub-eV range. 



6 Summary and Conclusions 



There is compelling evidence for neutrino flavor change driven by nonzero masses and mixing 
angles. Basically all oscillation data (with the only exception of LSND) are consistent within a 
three-neutrino framework. Within such framework, the global constraints from oscillation data 
can be summarized (see also Figs. 5 and 6 and Ref.^j) through the following ±2cr ranges (95% 
C.L.) for each parameter: 



Such ranges are marginalized over the four inequivalent cases [sign(±Am 2 )] <8> [cos<5 = ±1], i.e., 
over the two possible hierarchies and the two possible CP-conserving cases, which are currently 
undistinguishable. [Notice that the lower error on sin 2 #13 is purely formal, and corresponds to 
the positivity constraints sin 2 #13 > 0.] 

Concerning the observables sensitive to absolute masses (m«, turr and £), the situation is 
still unclear. Current constraints at the eV/sub-eV level are dominated by either upper bounds 
on £ from cosmology or by the 0^2/3 claim on mag, whose combination is not possible, however, 
at face value. Further studies and data are need to go beyond the general statement that 
neutrino masses should be smaller than ~ 1 eV, and to really explore the sub-eV range. 

Within the three-neutrino scenario, it appears that the most important unsolved problems 
require probing #13, S, the hierarchy, and the absolute neutrino masses. Needless to say, further 
experimental results or theoretical insights might also reserve big surprises and force us to go 
beyond such scenario, either by adding new neutrino states, or new interactions, or both. 
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Figure 1: Solar and KamLAND constraints in the mass-mixing plane (6m 2 , sin 2 #12) and for #13 

separately and in combination, at 1, 2, and 3cr level. 
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Figure 2: Three flavor analysis of solar and KamLAND data (both separately and in combination) in the parameter 
space (8m 2 , sin 2 612, sin 2 813). The contours represent projections of the region allowed at 2a (A^ 2 = 4). 



SK + K2K, ^,=0 



> 

0) 

7 
o 



£ 
< 



> 



E 

< 





1 


SK j 


'= 1 


.i.! 


1 = 



5 F- 

4 \ 

3 r 

2 i- 

1 I 



4 r 

5 \ 
2 I 
1 r 




K2K 




SK + K2K 



0.2 0.4 0.6 0.: 
sin 2 i!> 7 , 



Figure 3: Analysis of SK and K2K data (both separately and in combination) in the plane (Am 2 , sin 2 623) at 
813 = 0. The parameters (5m 2 , sin 2 812) have been fixed at their best-fit LMA values. 
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Figure 4: Three-neutrino analysis of SK+K2K+CHOOZ data, including subleading LMA effects. The results are 
shown as projections of the (Am 2 , s| 3 , s 2 3 ) allowed regions (at 1, 2, and 3<r), marginalized with respect to the 

four cases [cos 5 = ±1] eg) [sign(±Am 2 ) = ±1]. 




Figure 5: Global three-neutrino analysis of oscillation data. Bounds on s 2 3 are shown in terms of no = */ Ax 2 for 
KamLAND (dotted curve), solar (dot-dashed curve), solar+KamLAND (short-dashed curve), SK+K2K+CHOOZ 
(long-dashed curve) and all data combined (solid curve). In each case, the continuous parameters (Am 2 , S23, S13) 
and — if applicable — the discrete parameters [cos 8 = ±1] ® [sign(±Am 2 ) = ±1] are marginalized away. 




Figure 6: Global three-neutrino analysis of oscill ation data. Bounds on each of the parameters 5m 2 , Am 2 , 
sin 2 612, and sin 2 623 are shown in terms of no = U Ax 2 - In each plot, all parameters but the one in abscissa are 

marginalized away. 
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Figure 7: Analysis of oscillatory and non-oscillatory observables in the plane (?n/3/3 3 S). The 2<j horizontal band is 
preferred by the positive 0u2/3 claim, while the slanted 2a regions below (for normal hierarchy, N.H., and inverted 
hierarchy, I.H.) are preferred by all other data (i.e., oscillation plus cosmological data). The absence of overlap 
indicates tension among the data in the sub-eV range, as far as the standard three-neutrino framework is assumed. 



See also Rets.' 



